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The past few years have seen remarkable progress in the theory and phenomenology of QCD, 
bringing perturbative and nonperturbative methods into closer contact with each other and with 
experiment. 

Q\ '. I. INTRODUCTION: QCD A GROUP PORTRAIT 

G\ ■ 
OV 

t-H ' In this talk, I will summarize some recent developments in QCD, concentrating on, but not limited to, topics 
£>^' discussed at this meeting. Details, of course, can be found in the talks themselves, presented in lively sessions 
organized by Lance Dixon and Joey Huston for perturbative QCD, and by Paul Mackenzie and Claude Bernard 
^ \ for nonperturbative QCD. The future of the field is in the already-advanced convergence of these topics, sometimes 
thought of as nearly independent. Progress on heavy quark physics is described in the proceedings from separate 
parallel sessions and plenary talks. I must necessarily pass over some of the most interesting recent advances in other 
l , closely related fields as well for lack of space. 
^ | I have tried in the following to interleave perturbative and nonperturbative treatments of QCD dynamics. Let me 
begin with a few general comments on the place of QCD studies in high energy physics. 
1/-) ■ Why QCD? By now it is clear that QCD is a "correct", or phenomenologically relevant theory, at least the way 
classical Maxwell theory is "correct". Classical electromagnetism is an effective theory appropriate to the limit of 
many photons; quantum chromodynamics might be the long-distance limit of some more fundamental, underlying 
Q^ ■ theory. Its self-consistency, however, leaves us to free to study QCD in its own terms. 

This is a fascinating, and dauting, task, despite the fact that QCD is defined through a single dimensionless 
parameter, a s . The overall dimensional scale is determined by comparison with other interactions, for example, 
O || by measuring the strong coupling at the mass of the Z: a,(Mz). Its intrinsic interest aside, QCD has important 
"practical" applications, in the calculation of backgrounds to new physics and, for hadron colliders particularly, in 
predictions for new particle cross sections. From a theoretical point of view, however, what makes QCD so attractive 
• rH , is that it is a quantum field theory that requires all orders in perturbation theory and nonperturbative analysis to 
confront data at available energies. 

QCD is the exemplary quantum field theory. QCD exhibits most of the classic quantum field-theoretic phenomena 
discovered in the sixties and seventies, including asymptotic freedom, confinement, spontaneous symmetry breaking 
and instantons. The problems of strong interactions that gave rise to QCD were also, in the same time period, the 
original inspiration for concepts of duality and strings. As we saw at this conference, the strong interactions, now 
understood as QCD, are once more the meeting ground for field theory, duality and string theory. 

Again, QCD is evidently the correct theory of the strong interactions. Given its depth, however, "tests" of QCD 
should be thought of as tests, or perhaps better, "explorations" , of quantum field theory itself. Complementary to 
the extraordinary accuracy of selected perturbative predictions in quantum electrodynamics are the broad predictions 
of QCD, interweaving nonperturbative and perturbative scales and phenomena. 



*Plenary talk presented at the meeting of the American Physical Society Division of Particles and Fields (DPF 99), UCLA, 
Los Angeles, CA, 5-9 Jan 1999. 



1 



II. QCD AT THE SHORTEST DISTANCES 



It was the asymptotic freedom of QCD that first drew attention to gauge field theory as the unique description 
of the strong interactions at short distances. This theme continues to unfold in current experiments, and to provide 
a basis for extrapolations between energy scales and the detection of signals for new physics. Let us begin with a 
run-through of the underlying methods [j]. 

A. Methods 

Infrared safety. Infrared safe quantities are insensitive to long-distance effects, and may be calculated in perturbation 
theory An infrared safe quantity may or may not be directly observable. The classic examples include the total 

cross section for e + e _ annihilation to hadrons, and jet and event shape cross sections in e + e~. These can be written 
in the general form, 

Q 2 W(Q 2 ) = E c n(Q 2 /^ 2 ) «."(/*) . (!) 

n 

with the c n dimensionless functions of the ratio of the hard scale to the renormalization scale /j,. 

Factorization. Cross sections for deep-inelastic scattering (DIS) and for jet or heavy quark production in hadron- 
hadron scattering, are not purely perturbative, but appear as convolutions in parton momentum fractions of distri- 
butions / //i(CjM) of partons a in hadrons h, with perturbative hard-scattering functions <r PT , 

Q 2 o- phys (Q,x) = J2 °-a T (Q/^a s (n))®fa/h(n)= J2 I #^(s/£.Q/^««0*))/«/fc&P)> ( 2 ) 

partons a partons a 

where in the second equality we have exhibited the convolution appropriate to deeply inelastic scattering (with 
x = Q 2 /2p ■ q). Corrections to Eq. (||) are suppressed by 0(1/Q 2 ) 0. In this formula, there is a complementarity 
between the roles of parton distributions and the hard scattering. The distributions f a p t are universal among hard- 
scattering processes, but particular to hadron h, while the functions a PT are particular to the process, but universal 
among external hadrons. This last feature enables us to calculate realistic <r PT in "unrealistic" , but technically 
manageable (infrared regulated) scattering processes, in which the initial state hadrons are partons. 

Evolution. The physical cross sections of Eqs. ([!]) and @ above must both be independent of the scales \x that 
define the parton distributions: fj,da p hys/dfi = 0. This self-consistency requirement is readily translated into the 
"DGLAP" equation for the evolution of parton densities, 

an b 

Here, the dimensionless kernel P depends only on variables that are in common between the hard scattering functions 
and the parton distributions: a s and the momentum fractions. The scale-independence of physical quantities and 
their relations can be studied systematically ||. 

The idealized pattern for determining and applying the distributions may be summarized as follows. We measure 
one cross section, cr p h ys , a t momentum transfer Q . Given a "next-to- leading" order calculation of its hard scattering 
functions <ri NLO '', we determine NLO parton distributions fj^J^ ^ at /i = Qo- Using evolution, we can then predict 
Cphys for any hard process at all Q. 

The coefficients c n in Eq. (Q) are known for many processes to NLO ||. They have been determined at NNLO for 
inclusive DIS and Drcll-Yan B, and even to three loops S for selected quantities. Generally, two loop corrections are 



^1 have described some of the technical background in Ref 
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available only for single-scale processes, and the calculation of two-loop corrections for genuine scattering diagrams is 
an as-yet unsolved, but actively studied, problem in QCD HQ. 

Perturbative QCD is most successful for inclusive processes, and/or single-scale semi-inclusive. Evolution in DIS 
is perhaps still the best illustration (see below). The current experimental situation in hadronic hard-scattering is 
reviewed in Ref. In multiscale problems, logarithms of ratios of distinct but perturbative scales often require 
resummation to all orders. Formally beyond perturbative resummation, but not always less important numerically, 
are corrections suppressed by powers of the hard scale(s). In DIS, and a few other cases, these corrections can be 
described by the operator product expansion. We shall encounter below "generalizations" of this famous technique, 
usually in the form of effective field theories. 



B. Prime Examples 

Tevatron Jets. The most impressive success in orders of magnitude continues to be the Tevatron inclusive single-jet 
and dijet cross sections J| Jljjl , as illustrated in Fig. |[ which shows a plot from Ref. Jl^|. According to Eq. (||), these 
cross sections are of the form 

Vpp^J+X = fa/P ® a( alt°j+X ® fb/p ■ ( 4 ) 
ab 

We find a consistency between NLO theory and experiment at a few tens of percent, well within the overall systematic 
errors, over a range in which the cross section decreases by seven or so factors of ten. As the figure shows, reasonable 
choices of parton distributions (in this case CTEQ5HJ) can account even for the highest-momentum data, although 
a slight difference remains between the DO and CDF data at the high end. 
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FIG. 1. A comparison of single-jet inclusive cross sections to NLO theory from Ref. [[L2| 



DIS Scaling violations. Next, we should cite measurements of DIS structure functions [pli|-|l5[ for £ ± A — * i? ± A, 
through the cross sections 

' h > ' hV ° 2 [(l + (l-y) 2 )F 2 ( Xl Qi)~y?F L (x,Q 2 )T(l-(l-y) 2 )F 3 ( Xl Q 2 )} , (5) 



dxdQ 2 xQ 4 
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at HERA, the Tevatron, SLAC and elsewhere, with x — Q 2 /2p ■ q q 2 — —Q 2 y — Q 2 /xs. Because the evolution is 
universal in the factorized forms, 



F h (Q 2 ) = £ C a {Q 2 /n 2 ) ® f a/h (n 2 ) \ Q=/i 

a 

we may think of evolution as in, but not of, the nucleon, except perhaps in the smallest-x region, where target- 
dependent shadowing, in which partons begin to interfere with each other, comes into play. 

HI 97 preliminary 



• HI 97 prel. a HI 94-97 prel. QCDFit prel. 
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The quality of the data, and of QCD fits to it, is illustrated by recent results on F^cc, Q 2 ) from HI jl6), shown 
in Fig. H for various fixed x as a function of Q 2 . Notice in particular the sharp rise of F2 with Q 2 for the smaller 
values of x. This is a prediction of the evolution through parton branching described by Eq. (^). Deviations from 
such DGLAP evolution have been surprisingly difficult to find. The interesting results of a next-to-next to leading 
order analysis of DIS has been described in Ref. [|17j . At the other end of the spectrum, the excess of events reported 
two years ago at the largest x and Q 2 has all but disappeared JTsj ] . 

Multijets, and event shapes. In e + e _ annihilation, the Born level process e + e _ — > qq produces two jets, as is the 
case as well for hadron colliders, while in DIS, the Born scattering produces a single jet. In all three "canonical" 
scattering processes, the cross sections for the minimal reactions are known at NLO. Beyond this, NLO cross sections 
are available for two jets in DIS |L9]], for four in c + e _ ^0|, and NLO three-jet cross sections are just now becoming 
available for hadron- hadron scattering Recent years have seen an explosion of data on these processes, and 

new results were discussed at this conference in p|- |TT| , |Hi| , p2[ , along with jet fragmentation properties p^-|27j]. The 
production of vector bosons associated with jets [|28| is another important source of information about short-distance 
dynamics, and a reanalysis of W+l jet cross sections has brought them back into agreement with theory over the 
past year. 

Jets are theoretical-phenomenological constructs |2!| . Their value is not that they are an exact reflection of short- 
distance reactions, but that, if they are defined properly, they are related to them in a calculable fashion. Jets are 
generally defined in terms of energy flow into some angular region, or in terms of interative clustering schemes for 
the momenta of observed particles. Any (sufficiently smooth |50|) quantity that is insensitive to the emission of 
zero-momentum particles, or to the collinear branching of finite-energy massless particles, can be used to define an 
infrared-safe cross section. Event shapes in e + e~ annihilation or DIS are chosen for this property. The best-known is 
the thrust, defined for individual events as the maximum fractional projection of the momenta of observed particles 
on an axis, as that axis is varied about the unit sphere. Although event shape cross sections are infrared safe, they 
generally receive nonperturbative corrections that decrease rather slowly with energy, typically as a single power [ |3l| , 

a phys = <7 PT (l + 0(l/Q)) . (7) 

The theory of these power corrections is a passageway between the short-distance perturbative, and the long-distance 
nonperturbative aspects of QCD. We shall have more to say about them below. 

The a s lineup. Any of the short-distance cross sections above allows a determination of the strong coupling. For 
example, in a factorized cross section, Eq. (^|), given parton distributions /, we compare a FT (a s ) to experiment, 
and solve for a s , typically evaluated at a renormalization scale equal to the factorization scale. Other, more refined, 
choices || are related by perturbative corrections. For cross sections that are infrared safe, as jet or event shape cross 
sections at LEP 32|, the comparison is even more direct, although power corrections, as in Eq. (Q), must be taken into 



account for any precise determination |g^,g3[. Yet another way of determining a s is from lattice QCD calculations 
of energy level differences in heavy quark systems. These may be related, on the one hand to the strong coupling, 
and on the other to experiment, which then determines the size of the coupling at a scale that grows with the heavy 
quark mass. 
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For ease of comparison, the couplings can be evolved to the mass of the Z, at which a "world average" of a s (Mz) = 
0.119± 0.004 has been quoted in Ref. [Q. Since varied experiments measure the coupling over a wide range of scales, 
Fig. |^, which shows the variation of measured a s (/i) with /z illustrates one of the great successes of renormalized field 
theory. 



C. Resummation 

In the presence of multiple scales, high order contributions to hard scattering functions <r PT can become important. 
For example, in Eq. (0), the limit x/£ — > 1 is associated with integrable singularities that usually enhance the cross 
section. These "threshold" singularities were resummed long ago for the Drell-Yan cross section, and have been applied 
for some time to top production at leading- log accuracy |35|| . In QCD hard-scatterings, such as heavy quark and jet 
production, color exchange makes the resummation somewhat more complex beyond leading log. This problem has 
now been solved, and threshold resummation is understood in principle for a wide variety of hard scattering cross 
sections [p6|-^8|. As a practical matter, differing approaches to the inversion of certain Mellin transforms can lead to 
differing numerical predictions. From a broad perspective, however, the main lesson is that the theory stays relatively 
close to NLO for cross sections like top production, even in the presence of superficially large corrections at higher 
order. An important observation ]35| , |37| ], made particularly clear in the very recent publication |38| , is a marked 
decrease in factorization-scale dependence for resummed cross sections. 

Another example, currently being discussed widely, is related to the data of Fig. ||, single-photon and pion inclusive 
cross sections measured by the fixed-target Fermilab experiment E706 (3^^^]. So far, this data can be fit only by 
supplementing NLO with the old method of Ht smearing for the initial partons At the same time, it has been 
noted that the full range of fixed-target direct photon data may not be consistent among themselves ]4^ |. It should 
also be noted that at collider energies, experiment and NLO agree, at least at the higher transverse momenta [Q. 
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FIG. 4. E706 data versus NLO predictions and fcT-smearing. 



Can the resummation of higher orders in QCD lead to an effective kx smearing? For W or Z production at low 
transverse momentum, the answer is yes, and the formalism has existed for some time. In this case, fitting the cross 
section [ ff"H requires the introduction of nonperturbative parameters that are accessible to experiment. In the past 
year or two, there has been some discussion on the best way of going about this P5|,p6[ , but the underlying theory 
is relatively well-understood. The same cannot be said for W or photon cross sections at high transverse momenta, 
because the logarithms whose resummation requires the nonperturbative input at Qt ~ cancel order-by-order in 
perturbation theory in the calculation of a at high Qt- Interesting observations on the relationship between the two 
regimes have been made, however p7|, and futher progress in this direction can be anticipated in the coming year. 

One of the major challenges in perturbative QCD is the development of a theory of these and related higher-order 
effects, and a method for estimating their importance. We shall return to developments on resummation below. 



III. THE LONG AND SHORT OF HADRON STRUCTURE 
A. Parton Distributions 1999 

The parton distribution functions (PDFs) in Eq. (^|) summarize the structure of hadrons as seen at short distances, 
one parton at a time. As the hadron is probed at ever shorter distances, the distributions evolve perturbatively 
according to Eq. (|^), but the boundary condition for this evolution stands as a truely nonperturbative reflection of 
the dynamics that holds the hadron together. These parton distributions have been the subject of intense study since 
the recognition of approximate scaling in DIS structure functions thirty years ago. 

Over the past decade, two groups, CTEQ and MRS, whose memberships have themselves evolved somewhat, have 
undertaken coordinated "global" approaches to the determination of parton distributions, taking into account data 
from a variety of processes and momentum scales. The past year has seen the development of the latest, best fits of 
these two groups, MRST Q and CTEQ5 [Q, which were discussed and compared at the conference in flf9| . 

Global fits test the self-consistency of factorized cross sections in the sense that the fits are overconstrained, and 
because they can be checked against experiments not incorporated into the fits. Nevertheless, the fits must be improved 
as the data improves, and as it extends to extreme values of fractional momenta. Surprises can occur, especially in 
regions where a particular parton distribution does not contribute at leading order. Examples of such refinements 



7 



from 1998 involved the ratio of d to u quarks in the proton, as tested by the W asymmetry |50| and Drell-Yan fl5l| . 
Generally speaking, indirect constraints on parton distribution functions are tentative. 

Other cases where the results of global fits have been rethought involve higher-order or power ("higher-twist") 
corrections to the cross section. The extraction of PDFs requires a parameterization of such effects, assuming that 
they can be brought under control. Examples include the role of higher twist in the extraction of neutron PDFs from 
deuterium data p^ . 

Up to this cycle of global fits, the primary processes employed were DIS, Drell-Yan and direct photons, the latter 
thought to be especially valuable for constraining the gluon distribution, which does not appear at leading order in 
the other two. The data of E706 p9| , p0| , however, a sample of which was shown above in Fig. [|, has thrown this neat 
picture into disarray, as it disagrees decisively with the predictions of Eq. (||) at NLO. A kr smearing approach has 
been used in the MRST fits of last year p8| . The CTEQ5 fits abandon direct photons in favor of jet cross sections to 
constrain the gluon at large x [ ]T2[ |. 

Yet another interesting problem is the treatment of heavy quarks. For energies near a heavy-quark mass, toq, 
their production may be calculated directly as a hard process, part of <r PT in Eq. (^) . For energies much above toq , 
however, it may be advantageous to treat the heavy quark as a parton, thus automatically resumming logs of itiq/s. 
A number of schemes to make this transition have been proposed |l^,[l8|,^3],[34| , to treat charmed quark production at 
HERA. 

The sophistication of these considerations of global analysis, and the need to make accurate predictions at high 
energy raises the question of how to estimate uncertainties in the distributions |5f| . In part to explore these issues, 
new sets of distributions have been produced based on DIS data only |5(|, and methods have been introduced to 
quantify uncertainties systematically through statistical analysis [B7J. 



B. Spin and Off-diagonal Distributions 

The past few years has seen a rebirth of interest in the high-energy physics of spin, which, with advances 
in the technology of polarized beams and targets, has made possible the systematic study of spin at the par- 



ton level |13|. The DIS cross section for a nucleon with spin s may be represented in the standard form 



d a/dQdE — (a| M /2mQ )(E e / E'^L^W^, in terms of spin structure functions defined as 

= + 5 ' F/ ^AagV 91 (X, Q) + * , 6 M „A g g A [p ' ?*" - 8 • qV°\ 9*{x, Q) • (8) 

The function g\ has a particularly transparent interpretation in terms of quark helicity: <?i(a;, Q) = \ Ylf ejAqf(x, Q) + 
0(a s ), where Aqf is the difference between the distributions of quarks with helicity parallel to the hadron's helicity 
and against it. At this conference, precision data for g\ and g™ were presented by E155 [j58| . 
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One result of these measurements is a test of of the benchmark Bjorken sum rule, 



m(x)-9i(x))dx = - 
o 



9a 



9v 



(9) 



which is now verified to a new level of accuracy. 

The "spin" distributions of the nucleon, such as Ag^, do not necessarily describe its full spin content, and the 
possibility of orbital angular momentum must also be taken into account. At the same time, for gluons the distinction 
between these two types of angular momentum is not gauge invariant. This problem notwithstanding, an attractive 



formalism for the description of orbital angular momentum has been proposed ]^ 
5 L4 9 . g (0) + -Bg. g (0)], which arise in matrix elements of the energy-momentum tensor 



in terms of form factors, J g , 5 = 



(p + A| T% \p) = u(p + A) A q , g (A 2 ) -f? + B q , g (A 2 ) p" — A Q - (/i 



(10) 



with p — p + A/2. At zero momentum transfer, the J q g become expectation values of the angular momentum 
operators 



Jo = / d 3 x 



, J g = J d 3 x x x (E x B) . 



(11) 



The measurement of these off-diagonal matrix elements, unlike the diagonal matrix elements that define classic PDFs 
(see Eq. mw below), require the measurement of exclusive, or semi-exclusive amplitudes, such as "deeply-virtual 



Compton scattering", 7*(Q 2 ) + N(p) -> 7 + N(p'), with p 2 = p' 2 = 



'N- 



The measurement of such amplitudes 



is a challenge, but one that is of great interest for the Jefferson Laboratory facility, and there is a correspondingly 
vigorous theoretical program to study the factorization and evolution properties of off-diagonal distributions |30] [32| . 
Off-diagonal matrix elements interpolate, in some ways, between inclusive and exclusive processes, and between parton 
distributions and hadronic light-cone wave functions |63fl . 

The consideration of orbital angular momentum leads us to the doorstep of the long-distance, low-energy properties 
of hadrons, where progress has continued in lattice QCD and instanton studies. 







C. Lattice Hadron Spectra, Quark Masses 



Lattice methods |p4f approach QCD from a limit complementary to perturbation theory, and make possible direct 
calculation of long-distance properties of hadrons. Typically, lattice computations involve the evaluation of expecta- 
tions of nonlocal products of operators, such as 

Cj = (0 \J(x) J(0)| 0) , J = Hqi(x + A)O ijqj (x-A), (12) 

i 

with the qi quark fields, where Oij projects a set of quantum numbers. For x — ► oo, the x-dependence is dominated 
by the energy of the lowest-lying state(s) of the relevant quantum numbers. As noted above, such studies can be used 
to set the scale for a s by studying the hyperfine splitting in heavy quarkonia. 

The numerical evaluations of expectations like ( |l2| ) are said to be quenched, partially quenched, or fully unquenched, 
depending on how many species light fermions are allowed to fluctuate out of the vacuum. At the conference, progress 
was reported in quenched and unquenched lattice QCD, on the calculation of realistic spectra for hadrons, and of 
matrix elements of hadrons involving both heavy and light quarks p5| , |66| . Variations on this theme are now making 
possible the calculation of realistic decay matrix elements p6|-|68| . In addition, they allow the exploration of ideas 
on the mechanisms of confinement In alternative formulations, lattice methods may be applied to light-cone 

formulations of hadronic structure [ |63|]70|l . 

The current sophistication of lattice techniques, and the growing power of new machines, including Teraflop comput- 
ers at the U. of Tsukuba and at Brookhaven(RIKEN)-Columbia, are making possible a new generation of investigations 
of chiral symmetry breaking. 

Roughly speaking, chiral symmetry in QCD (also discussed at the conference in [[n]) reflects the observation that 
gluon emission doesn't change helicities. As a result, in the absence of quark masses in the propagator, left- and 
right-handed quarks decouple altogether in perturbation theory. This means that a massless quark stays massless 
in perturbation theory. Chiral symmetry, however, is broken nonperturbatively at zero temperature, but restored at 
finite termperatures, a transition which can also be studied on the lattice |^2|,^3|. At the same time, current algebra 
requires that the square of the pion mass would vanish linearly with the lightest quark masses, 

ml ~ mught ■ (13) 

Historically, this kind of relation has been difficult to realize on the lattice because of corrections that vanish as a low 
power of the lattice spacing in the continuum limit. 

With the advent of ever more powerful machines, the method of "domain wall fermions" [ [74[ makes possible new 
approaches to relations like Eq. (^). In the domain wall technique, the chiral symmetry is manifest, up to exponential 
corrections, at the price of introducing a 5th dimension, which we label by coordinate s, in the five-dimensional Dirac 
Lagrange density 

L = i> [j ■ D[A] + m(s)] V>, (14) 

where the mass parameter depends on s, and vanishes at endpoints and sq: m(0) = m(so) = 0. The light quarks 
are zero modes that propogate in "our world" at s = and sq, one wall for each helicity. These new methods have 



shown early success ]73| , |75j in preserving chiral symmetry, in terms of relations such as Eq. (13). In the lattice domain 
wall construction, the extra dimension is just a convenience, but it cannot be denied that the technique bears an eerie 
resemblance to the brane constructions of modern string theory p6[ . 



D. Hadrons and Instantons 

Successes in the treatment of chiral symmetry in lattice QCD lead us naturally to the reemergence of instanton 
studies of hadron dynamics M%. In the language of spontaneous symmetry breaking, the masses of hadrons are 
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related to the generation of quark "condensates" (qq) ~ + (qlQr) in the QCD vacuum, which couple (the 

perturbatively decoupled) left- and right-handed components of the quark field. 

Instantons may be thought of as tunneling events between the inequivalent QCD vacuum configurations that are 
distinuished by different phases of the nonabelian gauge fields at infinity (even at zero energy). The instantons 
couple to the quarks, and the process of tunneling produces an effective 2A^j-point interaction between their left- and 
right-handed components, of the general form 

(15) 

even though the quarks are massless. 

In the "instanton liquid" model, the nonperturbative soft gluon field is replaced by an ensemble of instantons and 
anti-instantons that couple left- and right-handed quark-antiquark pairs. Meson masses are produced as pairs "hop" 
between instantons, changing helicities as they go. In this model, mesons and baryons, with realistic spectra and 
matrix elements emerge, and provide a cross-check for lattice calculations j72[ at zero and at high temperatures. 
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IV. FACTORIZATIONS, EVOLUTIONS AND EFFECTIVE THEORIES 



Every hard-scattering experiment includes a complete evolution all the way from short distance to long distance 
dynamics. Factorization allows us to organize the long distance dynamics, and thus to calculate perturbative short- 
distance dependence, and compare the results to experiment. The essence of factorization is to interpret long-distance 
information in terms of matrix elements in the underlying theory. For example, in the classic factorization in Eq. (|^) 
for DIS structure functions of hadron h, F^(Q) = ^2 a C a (Q/ fj.) ® f a /h{l-i>), the quark (a = q) distributions may be 
interpreted as 

/,/*(*, m) = /_°° Z e ' ixPnv {h{p)l 5(^K7<i>nMMo) IMP)), (i6) 

where n is a lightlike vector not in the direction of the hadron momentum P, and /i enters as the renormalization scale 
for the matrix element, which is ultraviolet divergent for n 2 = 0. The function <& n (y,0) is a path-ordered exponential 
of the gauge field, 



®n(y,0) = -Pexp 



rv 

-ig I dX n ■ A(Xn) 



(17) 



which makes the matrix element gauge invariant. The evolution of Eq. ([?]) may be thought of as a consequence of 
the renormalization properties of the nonlocal operators in f a /h, which summarize an infinite set of twist- two matrix 
elements in the light-cone expansion. The same formalism is at the basis of heavy quark effective theory and of 
nonrelativistic QCD (NRQCD) ]7q |, although generally with a finite sum over operators rather than a convolution. 
In a sense, the operator $„ in Eq. ( Jl6| ) plays the role of the heavy quark field in heavy quark effective theory, as a 
nonrecoiling source of gluons. 

Over the past few years, factorization, analyzed in terms of effective operators, has been applied to multiscale hard 
scattering processes, with Sudakov |^,^0| and Regge high-energy limits. The former refers to cross sections with large 
momentum transfer and low QCD radiation, the latter to low momentum transfer and essentially unlimited radiation. 
The Regge limit is related to the total cross section, while the Sudakov limit highlights its short-distance components. 

One of the simplest examples of the Sudakov limit is the dijet cross section in e + e~ at fixed jet masses mf, i = 1, 2. 
In the limit of light jets, mf <C Q 2 , the dijet cross section is related to the integrated thrust cross section by 

1 , da{T) f 2 2 da ( m 2 + m 2 \ 

The cross section for "nearly-lightlike" jets in e + e~ satisfies a factorization J79| 
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da 



dmfdm 2 



Q 6 



f[Ji(p,Pi)®S(ji,pi), 



(19) 



up to corrections suppressed by powers of mf /Q 2 , with ft the 4- velocity of jet i. In (|T9|), there is a double factorization, 
separating the dynamics of the jets, included in the functions J;, from both the truely short-distance "coefficient" 
function C and from the dynamics of relatively low-energy partons emitted coherently by the jets and included in 
the function S. The "soft" function S is associated with a particularly interesting composite operator in QCD. S 
describes the emission of gluons whose wavelenths are so long that they cannot resolve the internal structure of the 
jets, and are thus generated by the product 



W(0) = $ A (00,0)^(00,0) 



(20) 



where, in the notation of Eq. fll7|), the <&'s are ordered exponentials, and ft is the velocity of the quark jet, and ft of 
the antiquark jet. 

We need not dwell on the nature of the convolutions denoted by ® in Eq. (|l9|), but the double factorization itself 



is adequate to imply a resummation 79 of double logarithms in 1 — T 



1 



ftot Jt 



, ,da(T') 



-2C F 



da 



l-T 



aQ 2 



die? a s (k 2 ^) 



(l ./ Q 2Q2 k T 



(21) 



where corrections include fewer logarithms of 1 — T in the exponent. These results are also related to the renormal- 
ization properties of the operators W in Eq. ( p(i| ) . 

Another application of factorization and effective operators is to resummation in the "Regge" limit s — > oo, t fixed, 
the "BFKL" regime for QCD. The BFKL equation may be derived from a "multiperipheral" reexpression of DIS 
factorization, Eq. (||) 



F(x,Q 2 ) 



C 



G(tQ 2 ) + 0{l/Q 2 ) 



= J d 2 k T c(^,Q,k T ^j^\k T ) +0(l/ln(l/aO) , 
where the A-r-dependent distribution i)j is related to the gluon PDF by p4|: 



G{i,Q 2 



d 2 k T ^,k T ). 



(22) 



(23) 



In the second form of Eq. (f22|), the roles of longitudinal and transverse momenta have been reversed, and corrections 
are suppressed only by logarithmics of x, rather than powers of Q. In these terms, the BFKL equation, 



dt 



d 2 k' T K,(k T ,k' T )^,k' T ) 



d 2 k T 



7T 2 J {k T -k' T f 



(24) 



with ip = (l/fcy describes the evolution of hr) in ^. The same equation may also be derived from the renor- 
malization of ordered exponentials, like Eq. ([l7]). Indeed, distributions of exponentials integrated over a transverse 
density pi 



J d 2 xx p(xt) P exp — ig J dx + A (x , xt) 



(25) 



are being studied to develop "unified" effective theories that describe the variety of evolution equations, including 
DGLAP, BFKL and others §1. 
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V. BFKL, DIFFRACTION AND COLOR DYNAMICS 



A. BFKL at NLO 

What's so special about the BFKL equation? It addresses the total cross section in gauge theory in terms of its 
fundamental quanta, the quarks and gluons. Now it is not entirely obvious that such a project will work, ultimately, 
in an asymptotically free theory with confinement, but if it does, it will say something fundamental about held theory. 
In an older language, this would be a theory of the "pomeron" ]8j|. Also, in the language of the parton model, 
BFKL appears to predict that, as we evolve to low x «-> high s 7 *jv at hxed Q, we reach a region of high parton 
density at nearly fixed (actually slowly diffusing [|S6||) virtuality. This is a new "intermediate" regime of QCD, between 
perturbative and hadronic phases p7[ . It is relevant to the screened "plasma" state, which we hope to encounter at 
RHIC. 

The solutions to the lowest order BFKL equation are of the form 



-ii/-l/2 



(26) 



The largest permissible value of ui gives the dominant low- a; behavior, which is found to be 



v = 0, to = ujq = 4AHn2 (a s /w) . 



(27) 



From the kinematic relation in DIS, s 7 *at = Q 2 (l — x)/x, the low-x behavior of ^{x,Ut) determines the large-s 
behavior of the 7*iV total cross section, for which the lowest-order BFKL result (^) gives 



' Ctot 



4AHn2 (a s /7r) 



(28) 



This is a derivation of Regge-like behavior for the total cross section from perturbative QCD, and, because it involves 
the exchange of no overall quantum numbers, may be considered as a perturbative model for the pomeron. 

Where should one look for BFKL behavior in experiment? Suggestions include correlations in dijet and rapidity-gap 
cross sections in DIS and pp, and in non-DGLAP evolution in DIS at low x and moderate Q 2 . In the first case, there 
may be hints in the dijet data from HERA and in the comparison of jet correlations at 630 and 1800 GeV jyj. In 
the later case, the strong rise in x of the structure functions Fi(x,Q 2 ) cannot be sustained indefinitely, since, by 
(P8|) , this would eventually violate unitarity bounds for the cross section. Before this happens, interference between 
partons, or "shadowing", which is absent in both DGLAP and BFKL evolutions, must begin to set in |pq| . One of the 
much-discussed data presentations of the past year, Fig. |^ from the ZEUS collaboration, shows the transition between 
perturbative and nonperturbative behavior in a particularly suggestive form |89[] . 

ZEUS 1995 Preliminary 



rvco O *l cn ^ c\ o ir 

= (vcv cv r> to — v 

o'o'o' o' & & & ~-'c\' 



ZEUS DATA 



DL 




10 10 ~ 10 10 ~ 10 10 



FIG. 6. Slopes of F2 as measured by ZEUS. 
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1998 was the year of the NLO BFKL kernel, the year in which the decade-long project of computing the next-to- 
leading order kernel in Eq. ( |24| ) bore fruit |xj. At NLO, K, is fairly complicated, but the effect on u> in Eq. ( p7| ) is 
simple enough, 



1 - 6.6 



(29) 



Now 6 and a half is not by itself a large number, but the size of this correction nevertheless presents a challenge, 
because unless a s is quite small, the second term may overwhelm the first and lead to an unrealistic falling cross 
section. In addition, it has been observed that the NLO kernel even implies a non-Regge behavior at high orders [ plj . 
Interesting responses to these challenges were discussed in |)2| . Evidently, the NLO fruit of BFKL will be an acquired 
taste, but the coming year surely promises intensive work and further clarification 



B. Diffraction and Diffractive PDFs 

Through the optical theorem, the total cross section is closely related to elastic and diffractive scattering amplitudes 
@. In diffractive DIS, we can relax inclusivity and thus probe QCD dynamics in the final state, while retaining a 
large momentum transfer. Convenient variables to describe diffraction are 

Ml + Q 2 x 

*-=wW /3 "' (30) 

where Mx is the mass of an observed system X moving in the "current" (i.e. photon) direction, and W 2 — s 7 * s . x-p 
is the fractional longitudinal momentum transfer from the nucleon to X, and [3 is the equivalent fractional momentum 
of a parton in the (hypothetical) exchanged "pomeron" . 

Diffractive events are typically defined by a large gap in rapidity between X and the elastically-scattered, or 
diffractively-excited, proton (or its low-mass fragments), which has experienced invariant momentum transfer t. In 
these terms, a fully differential cross section is 

In the simplest diffractive processes X consists of a single vector boson. In this case, the relevant amplitude is 
factorized in terms of off-diagonal PDF's, related for small x to the gluon distribution, which behaves as G(x) ~ x~ x . 
The resulting cross section is proportional to G 2 (x), and hence increases with W, as 

d(jD L U/4A(M X ) ( n 2] 

dM 2 x M x W ■ [62) 

There is compelling evidence for this behavior, with a value of A increasing with Mx, suggesting once again a Regge- 
like behavior reminiscent of BFKL. This connection has yet to be completely explored. 

High-Q 2 DIS diffractive cross sections may be factored using specifically diffractive PDFs, also referred to as 
fracture functions p4| , 

F D =Y, C a®fa- (33) 
a 

Phenomenological fits to f^„'s have been carried out fog}] , and lead to predictions whenever a factorization like Eq. 
( |33| ) applies |^6|. It is important to realize, however, that because diffractive PDFs are not fully inclusive, they depend 
on the details of evolution into the final state, in particular the likelyhood of the target proton staying together. This 
probability cannot be expected to be the same in pp cross sections, where the fragments of two initial-state hadrons 
pass through each other, as in DIS, where only a single hadron is involved. And indeed, studies J97j have shown that 
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diffraction at the Tevatron is much less likely than would be suggested by a direct generalization of Eq. (|33|) to this 
case with universal diffractive PDFs. 

Nevertheless, double-diffraction (double rapidity gap) jet production is seen at the Tevatron, and the jets show a 
standard parton-parton Er-dependence ]97|] , indicating that the short-distance process is independent of the long- 
distance evolution. This suggests that another factorization is possible in this case, and may shed light on diffractive 
dynamics. 



C. Color Dynamics 

Diffractive processes are naturally interpreted in terms of color-singlet exchange in the i-channel. Large momentum 
transfer processes, however, are described perturbatively in terms of single-gluon exchange, carrying octet quantum 
numbers. Over the past few years, there has been progress in understanding the relation between these two pictures. 

Although color is not observable, representations {l,q = 3,.g = 8...) are, at least in principle. In NRQCD, a 
factorization that includes the mixing of operators with differing color content has already led to valuable insights 
and phenomenological successes J78[]. For high-energy processes, rapidity gaps were long ago suggested by Bjorken as 
an ideal arena to study color exchange, with singlet exchange expected to produce an excess of events with very low 
interjet multiplicity. 

It is clear, however, that it is not possible to separate short- from long-distance color exchange uniquely, since gluons 
of all momenta carry the same color content, and models in which the color content of the final state is determined 
at the longest distances have had success |)8|. At the same time, energy flow into regions between two high-py jets 
is senstive to all time scales between 1/|pt| and 1/Aqcd- Energy flow at the shorter time scales is both perturbative 
and sensitive to the color content in the i-channel. This observation led |9^] to an analysis of dijet cross sections in 
terms of energy flow Q c into the interjet region, in the range Aqcd < ft € \f— t. The cross section at fixed Q c 
may be factorized, and logarithms of Q^/t resummed. This behavior is found from the renormalization properties of 
composite operators that are products of ordered exponentials, which generalize Eq. ( p0| ) to the 2^2 scattering of 
partons with color exchange (labelled here by f), 

4 2 
i— 3 i—1 

where T f is a matrix that couples the color of the incoming and outgoing ordered exponentials $> ni , representing the 
active partons of the hard scattering. These operators mix under renormalization and induce an evolution that tracks 
mixing in the color space as the scattering particles (including qq, qg, gg) evolve from short to long distances. This 
analysis offers a new set of predictions for px, energy and rapidity dependence of gap events, which can be tested at 
Run II of the Tevatron, and at the LHC. 



VI. POWER CORRECTIONS 

There has been considerable interest in power corrections to infrared safe quantities. As noted in Sec. 2 above, 
such power corrections are quite important in the phenomenology of jet cross sections and event shapes in e + e~ 
annihilation ]22] , |32] | . Behind this work is a hypothesis, that it is not necessary to model the details of hadronization 
to parameterize leading corrections to perturbation theory, and a hope, that plausible parameterizations inspired by 
perturbation theory will lead to useful insights at the perturbative-nonperturbative interface The hypothesis 

seems to be correct; whether the hope will be realized remains to be seen, but there are preliminary indications that 



it might be |100 | 



How does perturbation theory imply nonperturbation corrections? In the calculation of any IR safe quantity at 
NLO, we always encounter integrals of the general form 
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rQ 

I(a,p) = f(x) a s (Q 2 ) I dk k p , (35) 



with p > — 1, where k may be thought of as a gluon momentum scale, and f(x) represents the remaining (IR finite) 
dependence. Now in many cases (see Eq. ( |2l| ) for example) we can argue (or derive) that higher order corrections 
modify (H) to 

rQ 

I^ csum \a,p) = f(x) / dk a s (k 2 ) k p . . . . (36) 
Jq 

This result shows that the perturbative reexpansion in terms of a s (Q) is asymptotic, with high orders that grow as 
n\ at large n. This information is encoded in the singularity of the perturbative expression for a s (k 2 ) at k — Aqcd- 
We can further reinterpret this behavior by means of a Borcl transform, but the inverse transform will not be unique 
in any case. Taking a more practical approach, we retain the perturbative factor f(x) and the perturbative range in 
the k integral in Eq. (|36|), and simply replace the lower end of the integral, k < p\, for some fixed fix, by a parameter 



I 



mi 

dk a a {k 2 ) k v ^ a p {^{) , (37) 



Since the overall integral in Eq. rt39) behaves as Q p+1 this is automatically a power correction. Clearly the value of this 



approach depends on the assumption that a p {jjL\) is in some sense universal [100]. Evidently, this is true approximately. 



and this method has found applications in models for DIS higher twist [101]. As the notation suggests, the parameter 
a p is often thought of as a reflection of a universal, nonperturbative low-scale running coupling. This is suggested by 
( |36| ) above, where it is a moment of the lowest-order running a s (k 2 ). It has been argued that the relation is more 
general, and that higher orders of a(k 2 ) may be incorporated into a reconstructed effective coupling, defined through 
dispersion relations [ 102| , 103| ]. 



In interpreting these developments, it is important to keep in mind that the values of higher- twist parameters cannot 



be defined independently of perturbation theory [ 100 1 , and that they will change as new orders, or resummations, are 
computed. A striking example of this effect was illustrated by the NNLO analysis of |l7j], which reduced the size of 
higher twist contributions, relative to those found in fits based on NLO. 

Another interesting application of these ideas is to resummed event shapes, as in Eq. (pl|), where a replacement 
like Eq. (|37j ) leads to a simple shift |l04[| in the thrust (T = 1 — t) distribution, which vanishes as 1/Q, 



da PT (t) da PT (t - A/Q) , Q ( 1 ( ;is . 



dt dt \{tQ) 2 , 

with A a constant related to ao in (p71). In a somewhat more general approach, we may once again factorize soft gluon 



emission into the region between the two jets, and derive a convolution expression for the cross section [105 



daMt) • ^ de /(e) ^(^/Q) +0 (4 S ) 1 (39) 



dt J w dt \tQ 2 



where /(e), a "shape function", has a field-theoretic interpretation [105,106], which involves the composite operator 



of Eq. (po|). /(e) is Q-independent and summarizes all \/{tQ) n corrections, while Eq. ( |39| ) reduces to (p8h with the 



replacement /(e) = S(e — A). A fit [[105| to /(e) using the extensive thrust data at Q = Mz, and the perturbative 
resummation of Jsif] , faithfully predicts da/dT for a wide range of Q, as shown in Fig. [?]. Given the discussion of 
the foregoing section, by following this line of reasoning we may hope to relate event shapes in e + e~ annihilation 
to energy flow in hadronic hard-scattering cross sections. This relation remains unexplored, although it is certainly 



related to multiplicity and correlation studies of the final states in jet events [107 
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0.00 0.05 0.10 0.15 0.20 0.25 0.30 

FIG. 7. The comparison of thrust data with a fit at Q = Mz based on Eq. (^), at energies (from bottom to top): 
1 = 14, 22, 35, 44, 55, 91, 133, 161 GeV. 



VII. QCD AT HIGH TEMPERATURE AND BARYON NUMBER 

I have already mentioned the path from perturbative QCD to high parton density through BFKL evolution. Cer- 
tainly, these considerations are made more interesting by the pending turn-on of the RHIC accelerator at Brookhaven, 
at which nuclei will be collided at unprecented energies. This development has led to a fresh look at QCD in "extreme" 
conditions, long of relevance to studies of the early universe. 

Color Superconductivity. At high enough density and temperature, the long-distance interactions that lead to 
confinement in the normal, hadronic phase of QCD are screened, and in some sense nonsinglet degrees of freedom are 
freed. In fact, QCD is expected to have a possibly quite rich phase structure in the plane of temperature (T) and 
baryon density (B). 

An exciting exploration of these features of the theory is the fresh look at the long-standing conjecture of color 
superconductivity at large B and low T, in the light of the instanton liquid model referred to in Sec. HID. It was 



realized in Ref. [108] that the four-fermion effective Lagrange density in Eq. ( |15| ) produces an attractive potential 
between quarks, which can lead to a condensate of quark Cooper pairs at the Fermi surface, in a manner analogous 
to superconductors of electric current. In fact, the condensate can be driven by gluon exchange, but the energy gaps 



produced by instantons are much larger in most, but not all [ 109 1 , of parameter space. Although more likely to 



be relevant for neutron stars than nuclear collisions, these interesting considerations were clearly inspired by RHIC 



physics, which has led to an efflorescence of studies of the QCD B/T plane [ 11C , 111 [. 

Energy Loss. As a final example, I will very briefly refer to some recent considerations on a topic of direct interest 
to RHIC and hadron-nucleus scattering, energy loss of fast partons in dense media. High-energy partons travelling 
through matter (partons or hadrons) will scatter and radiate, and their evolution into the final state will be modified 



in some way. In sufficiently inclusive processes, these effects are (perhaps surprisingly) higher twist |112|. It is of 



interest, however, also to look at changes in radiation at transverse momentum scales set by the medium, rather 



than a hard scattering. Along these lines, recent work on the energy loss [113] experienced by a quark travelling 
through a medium over length (L) has identified a QCD analog of the famous Landau-Pomeranchuk effect in QED. 
This work analyzes "induced" gluon radiation at A <C &t <C (Q), where (Q) is the typical momentum transfer in a 
projectile-medium scattering. If the Debye length of the medium is short enough in a hot, dense medium, Q could be 
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perturbative. 

The amplitude for an emission q 
a diffusion in b, 



g at impact parameter b is denoted f(b, t). The effect of the scatterings is 



da 



f(b,t)~b (Q 2 ^Q2 ) exp 



(40) 



with known corrections to the Gaussion. For large (Q) , the cross section involves the interference between amplitudes 
where the radiation occurs at different positions along the path through the medium, leading to a power spectrum in 
path length and frequency given by 



dl 



dio dz 

This results in an energy loss per unit length that grows with length (! 



d 2 bf(b,t) f(b,0). 



dE a s N c L 
~dz A\~ 



da 
~dQ 2 



(41) 



(42) 



at mean free path A, which is a diagnostic for \Q 2 -^2j, potentially able to distinguish the composition of the medium, 
whether hadrons, Debye-screened plasma or something else. Here again, the transition from energy loss at relatively 
low transverse momentum, to more inclusive cross sections, should be an interesting one |114|. 



VIII. CONCLUSION 



The central conclusion of this little review is the variety and vitality of the work itself. Beyond this, the nature 
of our knowledge of QCD is such that important ideas and techniques only require the possibility that they can 
be tested to receive further theoretical development. The ongoing experiments at LEP, HERA and the Tevatron 
have already transformed perturbative QCD into a truely quantitative discipline. The pending RHIC accelerator has 
inspired creative theoretical developments. A strong QCD component to future high-energy projects is sure to be 
richly rewarded by insights into quantum field theory . 
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